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Abstract. A short introduction summarizes the experimental high energy field of neutrino oscilla- 
tions and the oscillation notations. Then the results obtained by experiments on atmospheric neutrino 
oscillations are summarized and discussed. At the beginning of the new century a number of long base- 
line neutrino beams became available. The results obtained by different long baseline neutrino oscillation 
experiments are considered, including the first possible event v ^ — » v T in appearance mode. Finally 
conclusions and perspectives are discussed. 
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1 Introduction 

A high energy primary cosmic ray (CR), proton or nucleus, interacts in the upper atmosphere producing 
a large number of charged pions and kaons, which decay yielding muons and z^'s; also the muons decay 
yielding neutrinos and v e . The ratio of the numbers of to v e is ~2 and N^/Njj — 1. These atmospheric 
neutrinos are produced at 10-20 km above ground, and they proceed towards the earth. 

Atmospheric neutrinos are suited for the study of neutrino oscillations, since they have energies from 
a fraction of GeV up to more than 100 GeV and they travel distances L from few tens of km up to 13000 
km; thus L/E v ranges from ~1 km/GeV to ~10 5 km/GeV. 

The early water Cherenkov detectors 1MB pQ and Kamiokande [2 reported anomalies in the ratio of 
muon to electron neutrinos, while the tracking calorimeters NUSEX [3j, Frejus [4 , and the Baksan [5] 
scintillator detector did not find any. In 1995 MACRO found a deficit for upthrougoing muons [B:. The 
Soudan 2 experiment [7 confirmed the ratio anomaly. 
^ ! In 1998 Soudan 2 [8 , MACRO [9] and SuperKamiokande (SK) [ID] reported deficits in the fluxes 

with respect to Monte Carlo (MC) predictions and distortions in the angular distributions; instead the 
v e distributions were in agreement with non oscillated MCs. These features were explained in terms of 
v v v T oscillations. 

The MC atmospheric neutrino flux was theorically computed by many authors in the mid 1990s [TT] 
and in the early 2000s [12 . The last computations had improvements, but also a new scale uncertainty. 

Several long baseline v beams became operational: KEK to Kamioka (250 km), NuMi from Fermilab 
to the Soudan mine (735 km) and CERN to Gran Sasso(GS) (CNGS) (730 km). The K2K and MINOS 
experiments obtained results in agreement with the atmospheric v results. In the CNGS beam the 
OPERA experiment reported recently a possible first candidate event o v T in appearance mode. 

Very recently the MINOS experiment reported a possible anomaly of 17^ with respect to v ^ and on a 
possible indication for a 4 th neutrino. 

Long baseline experiments are expected to give increasing contributions to neutrino physics. 
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2 Neutrino oscillations and masses 



If neutrinos have non-zero masses, one has to consider the 3 weak flavour eigenstates v e , v T and the 
3 mass eigenstates v\, v^, v$. The flavour eigenstates v\ are linear combinations of the mass eigenstates 
Vm (and viceversa) via the elements of the unitary mixing matrix Ui m : 

3 

Vl = ^ UimlSm (1) 
m= 1 

In the conventional parametrization the matrix U reads as follows: 

10 \ / ci3 s 13 e iS \ / ci2 512 \ 

C 23 523 1 -512 C 12 (2) 

-523 C 23 / V Sl3e- iS C 13 J \ 1/ 

with S12 = sin ^12, and similarly for the other sines and cosines. 

In the simple case of two flavour (y^ v T ) and two mass eigenstates (^2, ^3) one has 



l/fjt = h> 2 COS # 2 3 + ^3 Sin #23 
= -^2 Sin #23 + ^3 COS 9 23 



where #23 is the mixing angle. The probability for an initial to oscillate into a v T is 



(3) 



->■ ^) = sin 2 20 23 sin 2 ( ^-22- J (4) 

The survival probability of a beam is 

P{v, „„) ~ 1 - P(„ M -+ e/ t ) = 1 - sin 2 20 2 3 sin 2 j ^ 27 ^^ (5) 

where Am| 3 = m| — m| and L is the distance travelled by the neutrino from production to the detector. 

Most measurements are from disappearance experiments. The SNO results may be considered as 
coming from indirect appearance experiments. At higher energies the OPERA experiment is a direct 
appearance experiment, while other experiments may obtain indirect appearance signals. 



3 Atmospheric v experiments 

After the early indications of an anomaly in atmospheric z/s [U [2j [6], in 1998 Soudan 2, MACRO and 
SK provided strong indications in favour of v ^ <H> v T oscillations. Confirming results were presented by 
long baseline experiments. 

Soudan 2 used a modular fine grained tracking and showering calorimeter of 963 t located under- 
ground in the Soudan Gold mine in Minnesota. The detector was made of lm x lm x 2.5m modules, 
surrounded by an anticoincidence shield. The bulk of the mass was 1.6 mm thick corrugated steel sheets 
interleaved with drift tubes. The Soudan 2 double ratio for the zenith angle range -1 < cosQ < 1 was R' 
= (N i/ /N e )^ATA/(N /i /N e )MC = 0-68 ± O.llstat, consistent with v ^ oscillations [8 . 

MACRO detected upgoing v^s via CC interactions — ^ u; upgoing muons were identified with 
streamer tubes (for tracking) and scintillators (for time-of-flight measurements), see Fig. 1 left [9} [T3]. 
Events were classified as: Upthroughgoing muons from interactions in the rock below the detector of 
with (E) ~ 50 GeV. Data were compared with the predictions of the Bartol96 [II], FLUKA, HKKM01 
[T2] MCs, Fig. 2 left. The shapes of the angular distributions and the absolute values favoured v ^ 
oscillations with Am^ = 0.0025 eV 2 and maximal mixing. Low energy events. Semicontained upgoing 
muons (IU) come from interactions inside the lower apparatus. Up stopping muons (UGS) are due to 
external interactions yielding upgoing muons stopping in the detector; the semicontained downgoing 
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Figure 1: Left) Cross section of the MACRO detector and sketch of event topologies. Right) Schematic layout 
of the SuperKamiokande detector. 

muons (ID) are due to downgoing ^ M 's interacting in the lower detector [12 . An equal number of UGS 
and ID events is expected. The average parent neutrino energy is 2-3 GeV. The data are compared with 
MC predictions without oscillations and show a uniform deficit over the whole angular distribution with 
respect to the Bartol96 MC predictions. 

<H> v T vs v ^ <o> v sterile- Matter effects for v ^ with respect to sterile neutrinos, yield different 
total number and different zenith distributions of upthroughgoing muons. The measured ratio between 
the events with -1 < cosQ < 0.7 and with -0.4< cosO <0 was R me os= 1-38 compared to R r =1.61 and 
RsteriZe=2.03. Thus the v ^ Vsteriie oscillations are exluded at 99.8% C.L. [13] 

energy estimate by Multiple Coulomb Scattering (MCS) of upthroughgoing muons. The estimate 
was made through their MCS in the absorbers. The ratios Data/MC noosc as a function of (L/E^) are in 
agreement with v ^ <H> v T oscillations [T3] . 

In 2004 in order to reduce the effects of systematic uncertainties in the MCs, MACRO used the 
following three independent ratios (it was checked that all MCs yield the same predictions): 
(i) High Energy Data: zenith distribution ratio: Ri = N ver ^/N/ lor ; (ii) High Energy Data: energy 
measurement ratio: R2 = N/ ow /^high] (iii) Low Energy Data: R3 = (Data/MC)/c//(Data/MC)/ J D+[/G5- 
Fitting the 3 ratios to the v ^ <o> v T oscillation formulae, yielded sin 2 2#23= 1, Am^ 2.3 -10 -3 eV 2 . 
Using Bartol96, one adds the information on absolute fluxes: 

(iv) High energy data (systematic error ~17%): R4 = N meas /NMc- 

(v) Low energy semicontained muons (scale error 21%): R5 = N meas /NMC- 

The last 2 ratios leave the best fit values unchanged and improve the significance of the oscillations [13] . 

SuperKamiokande (SK) is a large cylindrical water Cherenkov detector of 39 m diameter and 41 
m height containing 50 kt of water (fiducial mass 22.5 kt); it is seen by 50-cm-diameter inner-facing 
phototubes (PMTs), Fig. 1 right. The 2 m thick outer layer of water acts as an anticoincidence; it is seen 
by smaller outward facing PMTs. The detector is located in the Kamioka mine, Japan. Atmospheric neu- 
trinos are detected in SK by measuring the Cherenkov light generated by the charged particles produced 
in the neutrino CC interactions with water nuclei. Thanks to the high PMT coverage, the experiment 
detects events of energies as low as ~5 MeV. 

The large detector mass allows to collect a high statistics sample of fully contained events (FC) up 
to ~5 GeV. The FC events yield rings of Cherenkov light on the PMTs. FC events are subdivided into 
sub- GeV and multi-GeV events, with energies below and above 1.33 GeV. FC events include only single- 
ring events, while multi-ring ones (MRING) are treated as a separate category. The partially contained 
events (PC) are CC interactions with vertex within the fiducial volume and at least a charged particle, 
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Fi gure 2: Left) MACRO upthroughgoing muons compared with the oscillated MC predictions of Bartol96 (solid 
curve), HKKM01 (dash-dotted line), FLUKA fitted to the new CR fit (dashed curve) and FLUKA with the old 
CR fit (dotted curve). Right) ID events (black points) compared with the non oscillated Bartol96 MC (dashed 
line) and <H> v T predictions. 

typically the /i, exits the detector (the light pattern is a filled circle). Upward- going muons (UPMU), 
produced by v ^ interacting in the rock before the detector, are subdivided into stopping ((E) ~7 GeV) 
and throughgoing muons ((E) ~ 70 80 GeV). 

Particle identification is performed using likelihood functions to parametrize the sharpness of the 
Cherenkov rings, which are more diffused for electrons than for muons. The zenith angle distributions 
for e-like and /x-like sub- GeV and multi-GeV events are shown in Fig. 3 left. A MC problem exists also 
in SK: to reduce these problems the normalization is now left as a free parameter. 

The ratios e-like events/MC do not depend from L/E^ while /i-like events/MC show a dependence 
on L/Ej, consistent with the oscillation hypothesis, Fig. 3 right. In 1998 the overall best fit of SK data 
yielded maximal mixing and Am^ 3 = 2.3 • 10 3 eV 2 (2.5 in 2004). In 2010 SK started fits with three 
neutrino types: they obtained maximal mixing for #23 and Am^ ^2.2-2.3 • 10 3 eV 2 [T4] . 

Exotic oscillations. MACRO and SuperK data were used to search for sub-dominant oscillations 
due to Lorentz invariance violation (or violation of the equivalence principle). In the first case there could 
be mixing between flavor and velocity eigenstates. 90% C.L. limits were placed in the Lorentz invariance 
violation parameters | Av | <6 • 10~ 24 at sin 2 2<9 v =0 and | Av | <4 • 10~ 26 at sin 2 2(9 v =±l [T5l H5| fTT] . 

Neutrino decay could be another exotic explanation for neutrino disappearence; SK data do not 
favour v decay and no radiative decay was observed in specific experiments [18]. 

New limits were obtained on proton decay and on GUT magnetic monopoles by SK and MACRO. 

4 Long base line v experiments 

Atmospheric neutrino oscillations have opened a new window into phenomena beyond the Standard Model 
of particle physics. Long baseline neutrino experiments allow further insight into v physics. The first 
long baseline v beam was the KEK to Kamioka beam (250 km), the 2 nd was the Fermilab to the Soudan 
mine beam (NuMi) (735 km) and the third was the CERN to Gran Sasso (CNGS) (730 km). 
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Figure 3: Left) SK zenith distributions (black points) for e-like and /x-like sub-GeV, multi-GeV, throughgoing 
and stopping muons. Solid lines are no oscillation MC predictions. Right) L/E^ distribution for /x-like events. 



The K2K experiment, which used SuperKamiokande as the far detector and the low energy beam 
from KEK to Kamioka, confirmed the results of atmospheric v ^ oscillations [19 . 

The MINOS experiment on the NuMi low energy neutrino beam is a large magnetised steel scintillator 
tracking calorimeter, complemented by a similar near detector and a calibration detector, Fig 4 left. 
The experiment obtained important results which confirmed the atmospheric v oscillation picture with 
maximal mixing and Am^ 2.35 eV 2 , Fig. 4 right [20]. Recently MINOS reported a possible 40% 
difference between and oscillation parameters, Fig 4 right [21 . They plan to check this point with 
new runs, in particular with the z7^ beam. The MINOS experiment obtained also several interesting limits 
on sterile neutrinos, on Lorentz invariance and measured the neutrino velocity and the ratio for 
atmospheric muons [22]. The same ratio was also measured by the OPERA experiment [23], Fig 7. 

The CNGS v ^ beam. [23]. A 400 GeV beam is extracted from the SPS and is transported to the 
target. Secondary pions and kaons are focused into a parallel beam by 2 magnetic lenses, called horn and 
reflector. Pions and kaons decay into v ^ and \i in a long decay pipe. The /i's are monitored in 2 detectors. 
The mean energy is 17 GeV, the 77^ contamination ~2%, the v e (V^) <1% and the number of v T is 
negligible. The muon beam size at the 2 nd muon detector at CERN has a ~1 m; this gives a beam size 
at GS of (7 ~1 km. The first low intensity test beam was sent to GS in 2006 and 3 detectors (OPERA, 
LVD and Borexino) obtained their first events. The shared SPS beam sent a pulse of 2 neutrino bursts, 
each of 10.5 /isec duration, separated by 50 ms, every 12 s. The v ^ beam operated regularly in 2008, 
2009 and 2010. 

LVD is an array of liquid scintillators with a total mass of 1000 t, designed to search and study T^T's 
from gravitational stellar collapses. It monitors also the CNGS neutrino beam [25] . 

The ICARUS T600 detector is an innovative TPC liquid argon detector installed in Hall B of 
Gran Sasso [24]. It represents an intermediate technical step toward a more massive detector and it can 
offer several interesting possibilities. Now it is studying neutrino interactions in the CNGS beam with a 
fiducial mass of 480 t [26] . 

OPERA, in Hall C at Gran Sasso, is a hybrid emulsion-electronic detector, designed to search for 
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Figure 4: Left: photograph of the MINOS far detector: A indicates the magnetized iron sheets, C is the electric 
coil, B is an anti coincidence scintillator system. Right: MINOS determination of Am| 3 and sin 2 2^23 for and 
7^7 oscillations. 

^ —> v T oscillations in appearance mode [27]. The v T appearance is made by direct detection of the 
r lepton from v T CC interactions and the r lepton decay products. To observe the decays, a spatial 
resolution of 1 jam is necessary; this is obtained in emulsion sheets interleaved with thin lead target 
plates (Emulsion Cloud Chamber (ECC)). The OPERA detector, Figure 5, is made of two identical 
supermodules, each consisting of a target section with 31 target planes of lead/emulsion modules (bricks), 
of a scintillator tracker detector and a muon spectrometer. 




Figure 5: Schematic of the OPERA experiment at Gran Sasso 



The muon spectrometer consists of 2 iron magnets instrumented with Resistive Plate Chambers (RPC) 
and precision drift tubes. Each magnet is an 8 x 8 m 2 dipole with a field of 1.52 T in the upward direction 
on one side and in the downward direction on the other side. Precision drift tubes, measure the muon 
track coordinates in the horizontal plane. The muon spectrometer has a Ap/p^0.25 for muon momenta 
of ~25 GeV/c. The basic target module is a brick consisting of 56 lead plates (1 mm thick) and 57 
emulsion layers. A brick has a size of 10.2 x 12.7 cm 2 , a depth of 7.5 cm and a weight of 8.3 kg. Two 
additional emulsion sheets, the changeable sheets (CS), are glued on its downstream face. Within a brick, 
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Figure 6: Display of the first — > v T candidate event from the OPERA experiment. 



the achieved space resolution is <1 /im and the angular resolution is ~2 mrad. Walls of Target Tracker 
Scintillators provide the v interaction trigger and identify the brick in which the interaction took place. 
The bricks were made by the Brick Assembly Machine (BAM) and are handled by the Brick Manipulator 
System (BMS). A fast automated scanning system with a scanning speed of ~20 cm 2 /h per emulsion 
(each 44 /am thick) is needed to cope with the analyses of so many emulsions. For this purpose were 
developed the European Scanning System [28] and the Japanese S UTS [29] . An emulsion is placed on a 
holder and the flatness is guaranteed by a vacuum system. By adjusting the focal plane of the objective, 
16 tomographic images of each field of view are taken at equally spaced depths. The images are digitized, 
sent to a vision processor, analyzed to recognize sequences of aligned grains. The 3-dimensional structure 
of a track in an emulsion layer (microtrack) is reconstructed by combining clusters belonging to images 
at different levels. Each microtrack pair is connected across a plastic base to form the basetrack. A set 
of base tracks forms a volume track. 

Recently OPERA found a first candidate event — ^ v T [30] . We recall the main phases of the 
measurements. The electronic detector found one event which had the features of a NC interaction, 
that is without a muon. The Target Tracker scintillators individuated the brick in which the interaction 
occurred. The CS of the brick were removed and developed in the underground lab. Their analysis 
corroborated the TT findings and gave the precise location of some tracks. The brick was then brought 
to the outside lab, exposed to cosmic ray muons (in order to have more precise emulsion alignments) 
and developed. The emulsion sheets were later shipped to one of the outside laboratories of the col- 
laboration. There, measurements were made according to the following scheme: starting from the CS, 
tracks were reconstructed in the "scan back" procedure. The vertex of the interaction was located; the 
measurements proceeded with the "total scan" in one cm 2 around the vertex, which found several tracks. 
This was followed by the "decay search" and "scan forth" procedures to better determine the track with 
a kink and measure all the tracks of the event. Other bricks were extracted and measured to check for 
interactions, for black tracks of nuclei and for gamma rays. The event was visually checked and was 
sent to a second lab for a repetition of the measurements. The averages of the 2 measurements were 
used for track kinematics analysis, determination of the momentum of each track by multiple Coulomb 
scattering methods, and kinematic analysis of the candidate event, Fig. 6. The decay of the short track 
was compatible with r~ —> p~v T . The background from many sources was estimated. The probability 
that the background fluctuates to one event is few percent giving a statistical significance of slightly over 
2<j for the observation of 1 real event. One now waits for further events! 
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Figure 7: Compilation of ratios = /N^- in cosmic rays measured at ground level and underground. 
Notice in particular the data obtained at large e^cos^* by MINOS, LVD and OPERA. 



5 Conclusions. Outlook 

The atmospheric neutrino anomaly became in 1998 the atmospheric neutrino oscillation hypothesis, which 
was later confirmed with more data and from the new data from the first two long baseline neutrino 
experiments. All the experiments were disappearance experiments. All agree with maximal mixing and 
with the Am^ values indicated below in the 2 neutrino oscillation scenario (for SK also in the 3 neutrino 
scenario [14 ): 

- In 1998: 2.3 (SK), 2.5 (MACRO) 

- In 2004: 2.5 (SK), 2.3 (MACRO) 

- In 2006: 2.8 (K2K) 

- In 2008: 2.43±0.12 (MINOS) 

- In 2010: 2.35±g;£g (MINOS),.. .(2.2-2.3)(SK 3-flavors)(2.51 and 2.11 at the Neutrino 2010 Conference). 
The variations in time and the differences among the values from the main experiments are small and 
well within their experimental uncertainties. 

Now there is also one event from the OPERA appearance experiment and one is waiting for more 
such events. 

Solar and atmospheric neutrino oscillations provided the first indication for physics beyond the Stan- 
dard Model of Particle Physics [31 . One needs further data to clear several points, like the recent results 
of MINOS and the results from the experiments MiniBooNE [32 and LSND [33]. 

New experiments, like T2K [3l] and Nova [35], should improve considerably the situation, also using 
new particle production data [36]. The main aim of the new experiments is the determination of the 
mixing angle #13 and the search for a possible CP violation in the leptonic sector [37J [38] . 

Acknowledgments. I acknowledge many colleagues at CERN, LNGS and Bologna for discussions 
and advices. I thank drs. M. Errico and R. Giacomelli for technical help. 
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